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Abstract 

The Drude weight D and the dc-conductivity <r^ c (T) of strongly correlated 
electrons are investigated theoretically. Analytic results are derived for the 
homogeneous phase of the Hubbard model in d = oo dimensions, and for 
spinless fermions in this limit with l/d-corrections systematically included to 
Ph lowest order. It is found that <r^ c (T) is finite for all T > 0, displaying Fermi 

liquid behavior, <r^ c oc 1/T 2 , at low temperatures. The validity of this result 
for finite dimensions is examined by investigating the importance of Umklapp 
[- — scattering processes and vertex corrections. A finite dc-conductivity for T > 

^ is argued to be a generic feature of correlated lattice electrons in not too low 

dimensions. 
^ PACS: 71.27 +a, 72.10-d 
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I. INTRODUCTION 



Strongly correlated electron systems, such as heavy fermions and high-T c superconduc- 
tors, have very unusual thermodynamic and transport properties. To understand the physi- 
cal mechanisms underlying these properties, the temperature-dependent dynamic ("optical") 
conductivity o~(u, T) is one of the most meaningful physical quantities to be studied both 
experimentally and theoretically. Of special interest is the low-frequency behavior of the 
conductivity since it provides information about scattering processes close to the Fermi 
level. Theoretical studies of a(ui, T) for u — >• so far concentrated on the limit T = [1-8] 
where the real part of the conductivity can be written as 

lini Re a(cu,T) = D8(u) + <T incoll (a;) • (1) 

Here D is the weight of the Drude peak and is an incoherent background. As pointed 

out by Kohn [1] the Drude weight D is of particular importance since it provides a direct 
and sensitive criterion for a metal-insulator transition. 
The behavior of the dc-conductivity 

a dc (T) = a(u = 0,T) (2) 

at T > is not so well understood theoretically and received attention only most recently 
[9-13]. If a system has continuous translational invariance a^c must be infinite even at 
finite temperatures since the momentum acquired from the external electric field cannot 
be degraded (here we do not consider scattering with phonons or impurities). Indeed, the 
Hamiltonian H commutes with the total momentum P operator 

[H, P] = 0. (3) 

Here and in the following operators carry a hat. Hence any expectation value of P and thus 
of the current operator J = — eP/(mLi), is constant in time (here the electron charge is 
— e, m is the electron mass and L 1 is the length of the sample in the direction of P). No 
damping occurs. 

For a lattice system, characterized by a tight-binding Hamiltonian, this argument no 
longer holds. The current operator and the Hamiltonian do not commute and there is no 
reason why the conductivity should not be finite. Yet even on a lattice there are constraints 
on the destruction of momentum: the crystal momentum is conserved modulo reciprocal 
lattice vectors. In the absence of phonons or impurities crystal momentum can only be 
destroyed by Umklapp scattering. Now, it is not a priori clear whether Umklapp processes 
alone are sufficient to reduce a dc (T) to a finite value. For one-dimensional fermion systems 
this problem has been investigated in detail by Giamarchi [9]. In this paper we will approach 
the question of whether or not the dc-conductivity of correlated electrons may diverge at 
T > from the opposite direction, i.e. the limit of high dimensions. 
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In the following we consider a general lattice Hamiltonian of the form 

H = E c JCT + I E (4) 



^11 7,7 *J 



where cf a (c^) is the creation (annihilation) operator of a electron with spin a at site i 
and hi = ^ta^ia i s the total particle density at site i. The coefficients Uj and C/^ are 
hopping and interaction matrix elements, respectively. The coordination numbers Z^ count 
the number of sites which are in the same point-symmetry group as j with respect to site i. 
We assume the system to be invariant under discrete lattice translations and transformations 
of the point-symmetry groups of each site. In all explicit calculations we restrict ourselves to 
nearest-neighbor hopping on a hypercubic lattice (with coordination number Z = 2d, where 
d is the dimension), i.e. assume Uj = —t* for nearest-neighbors and Uj = otherwise. We 
set t* = 1, i.e. measure all energies in units of t*. The kinetic energy then reads 

Hkin = — L E = E e kHo ( 5a ) 

V Z ft)," k.a 



where er is the energy dispersion, 



er = = E l * os k< ( 5b ) 



Vzu 



with unit lattice spacing, and h^ denotes the moment distribution operator. Due to 
particle-hole symmetry it is sufficient to discuss densities n = N/L not larger than half- 
filling. Here iV is the total number of particles and L is the number of lattice sites. The 
doping 8 is measured relative to half-filling, i.e. 8 = 1 — n for the Hubbard model and 
8 = \ — n for spinless fermions. 

The scaling of the hopping by the square root of the coordination numbers in (4) is 
essential to obtain a well-defined limit d^oo for fermionic lattice models [14]. The scaling 
of the interaction terms in (4) is more straightforward: it must be chosen such that the 
total interaction strength over all sites does not diverge [15]. For explicit calculations we 
will consider electrons (with spin a =f , J.) interacting via a Hubbard interaction 



HHubbard 



-ni E + U^h^hn, (6a) 



as well as spinless fermions (no spin index) with nearest-neighbor interaction 

H S pi n l ess = -j= E <^~Cj + ^ fliflj . (6b) 

* (ij) (ij) 

The limit d — > oo is carried out on the level of ordinary diagrammatic perturbation theory 
where it leads to considerable simplifications [16]. These simplifications make an analytic so- 
lution of (6b) possible since the Hartree approximation becomes exact in this case [15,17-19]. 
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By contrast, the Hubbard model remains dynamic even in d = oo [15] and thus cannot be 
solved analytically. It reduces to a single-site theory [22] which is equivalent to an effective 
single-impurity problem with a self-consistency condition [23,24]. This equivalence makes 
numerical [24-27] and approximate [23,27] solutions possible. 

In the following we will investigate the dc-conductivity a dc (T) for the two models (6a,b) 
in the limit of large dimensions d, i.e. large coordination number Z (this limit is denoted 
by d — > oo or, equivalently, by Z — ► oo). In the case of the Hubbard model we work 
strictly at Z = oo. However, in the case of spinless fermions all terms up to order 1/Z are 
explicitly included. Not only is this managable in this case [19], but - more importantly - it 
is mandatory for obtaining finite results for a: in the limit Z = oo the theory is purely static 
(Hartree theory), leading to an infinite conductivity, since scattering processes are absent. 
The latter first enter into the theory through the inclusion of 1/Z-corrections [28], whereby 
genuine correlation terms are generated. We will refer to this limit as "spinless fermions 
with 1/Z-terms". 

To be able to study the effects of correlations on cr rfc (T) at low temperatures we will 
deliberately ignore the onset of long-range order (i.e. spin and charge order in the Hubbard 
model and spinless fermion model, respectively), since the existence of an energy gap would 
naturally lead to an exponential suppression of a dc (T) for T — > 0. 

The paper is structured as follows: in Section II the dynamic conductivity a(co,T) is 
calculated in the limit of high dimensions. A detailed investigation of the limit oo — ► 0, i.e. 
of the Drude weight D at T = and the dc-conductivity a dc (T) at T > 0, is presented 
in Section III. Explicit results are obtained for the Hubbard model in d = oo and spinless 
fermions with 1/Z-terms. The finiteness of a dc (T) at T > and its validity for dimensions 
d < oo is discussed in Section IV by examining the importance of Umklapp scattering 
processes and vertex corrections. 

II. CONDUCTIVITY a(v,T) IN THE LARGE- £> LIMIT 

The dynamic conductivity can be calculated from the Kubo formula as 

a(u,T) =g[a 1 {u),T) + a 2 {u),T)\ (7a) 

where 

, _ x ie 2 r d d k d 2 £r . , , 

"W^/m^ 1 (7b) 

BZ v ' 

a 2 (co,T) = — X j \co,T). (7c) 

CO 

Here g is a spin degeneracy factor, i.e. g = 1 for spinless fermions and g = 2 for the 
Hubbard model and x JJ (o;) is the retarded current-current correlation function; the current 
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is chosen to flow in the 1-direction. Since the current flows into a given direction, i.e. is 
proportional to the effective hopping matrix-element which itself is of order \j\[Z small, the 
current-current correlation function, and thereby the conductivity a itself, is proportional 
to the square of the effective hopping. Hence it is of order l/Z small. Strictly speaking the 
scaling of the hopping would therefore imply that the conductivity - like any other transport 
coefficient - would vanish in the limit Z — > oo. In the following we therefore investigate the 
leading order contribution in l/Z to a, i.e. lim^oo Za(u,T) which essentially represents 
the sum of that quantity over an entire shell of nearest-neighbor sites. 

For the infinite dimensional Hubbard model the conductivity is given diagrammatically 
by the dressed bubble-diagram [29,30,10,18]. Vertex corrections do not enter since they 
are of order l/Z 2 . The conductivity is hence given by a convolution of two one-particle 
propagators G^(ico n ) = G^ n = [iu n + ji — S m — e^] -1 where S m = T,(ico n ) is the self-energy 
(which is ^-independent in d = oo [14,15]) and /i is the chemical potential, i.e. 



AT r d^k 

X n (iu n ,T) = -y £/ 7^ sin2 (fci)%, m G *A«M»-", • ( 8a ) 

m,Z BZ 



For spinless fermions with 1/Z-terms the r.h.s. of (8a) yields the auxiliary function 
xi\ioi) n ,T) from which the full correlation function follows as [28] 

In the limit d — >• oo the term sm 2 (k 1 ) can be replaced by its average in the Brillouin zone, 
i.e. by 1/2. Then the integrand in (8a) depends on k only through the energy dispersion 
of the non-interacting electrons. This demonstrates clearly the mean-field character of the 
d = oo limit. The density of states (DOS) of the non-interacting system, N (e), which for a 
hypercubic lattice (see (5b)) in d = oo takes the form [14] 

1 exp(- £ 2 /2), (9) 



27T 

may then be used to express (8a) as an energy integral 

2T 

m,l 



2T r 

(iu n ,T) = J deN (E)G m (E)G l (e)6^ m . Wl , (10) 



-oo 



where G n (e) = [ico n + ji — T, n — e] 1 . If the self-energy S m is known the best way to proceed 
is to use a partial fraction expansion 

G^G^e) = -[{iuj m - S m ) - (m - S / )]- 1 [G™(£) - G,(e)] (11) 

which permits one to perform the ^-integral in (10) as 
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where g m = j d£N {e)G m {e) is the local one-particle propagator. Analytic continuation then 
leads to 

x n (uj,T) = ^(B 1 + B 2 ), (13) 

where [28] 

— oo 

V ' _^ M ; \w-[S R (w' + a;)-£ R (a/)] w - [S A (W) - S A (W - w)] J V j 

Here g R and g A are the retarded and advanced local one-particle propagators, respectively, 
and f(x) = (e@ x + 1) _1 is the Fermi function. In the limit T, to — > the term B\ dominates 
since it reduces to a constant and hence leads to a divergent contribution to <t 2 (w, T), (7c), 
and thereby to a(u,T) itself. However, in the case of spinless fermions with 1/Z-terms B 2 
has to be taken into account, too (see below). 



III. THE LIMIT to 

A. T = 0: study of the Drude weight D 

To obtain an analytic expression for the Drude weight D in (1) we use (14a) for B 1 in 
the limit of small to 

BAu, 0) = 2mN(0)- . (15) 

v ' v ' (m*/m)LO + 17] v ' 

Here iV(0), the value of the DOS N(u) = — / K~ 1 lmg R of the interacting electrons at the 
Fermi level, enters. Furthermore m* , with m* jm = 1 — dReY</duj\ w = , is the effective mass, 
and it] is an imaginary part which is infinitesimally small at T = but becomes finite at 
finite temperatures. From (7a) we find for D 

D = 2 T dcu Recr 2 (w,0) (16) 

J —rj' 

where rj' is also infinitesimally small at T = 0. The contribution of 01, (7b), need not be 
considered here since it is purely imaginary. For the same reason B 2 is unimportant in the 
evaluation of D for the Hubbard model. Combination of (15) and (16) yields 

47re 2 

ZD Hubbard = — — N(0). (17) 
m* jm 
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The situation is different in the case of spinless fermions with 1/Z-terms [28]. Here the 
limiting value of B 2 = Trfy -1 (Hkin) for small ui must be taken into account, too, where 

7 = 1- y/ZZ Fock (18) 

is a renormalization factor resulting from the frequency-independent Fock-contribution to 
the self energy [31]. Furthermore, N(u) must be replaced by the spectral density 

N c (lo) = --Im c(lo) (19a) 

defined through 

c(u) = r de K i £ ) - (19b) 

with N c $(e) determined from N (uj) as 

^7 = 2 ^ (19c) 

In d = oo, with N (ui) given by (9), this implies N^° (e) = ^N^°(e) and therefore 
cr,a{u) = 29 RA^) an d N c (e) = 2^( £ )' Finally? by increasing the effective hopping 
of spinless fermions with 1/Z-terms [31] the factor 7 also reduces the effective mass, i.e. 
m* /m = (1 — dReT,/dco\ UJ ^o)/'j. In this way one finds in analogy to (17) 

7D 4reVV c (0) , 2m 

spmZ ^ m*/m + 2Z- 1 [/7V c (0)' 1 J 

For the evaluation of (17) and (20) at T = it is helpful to note that in d = 00 the quantities 
/i — S(0) (for the Hubbard model) and [/i — S(0)]/7 (for spinless fermions with 1/Z-terms) 
are independent of U. This is a direct consequence of Luttinger's theorem on the constancy 
of the volume within the Fermi surface [32] and the fact that S is strictly local (Hubbard 
model). It even holds when S contains terms which are proportional to the kinetic energy 
(as in the case of T 1 Fork for spinless fermions with 1/Z-terms), even when they are frequency 
dependent as in the case of the Hubbard model with 1/Z corrections. Consequently, N(0) 
in (17) and 7iV c (0) in (20) may be replaced by iV (0) and 7V Ci0 (0), respectively, i.e. by the 
quantities at U = 0, at the same filling. 

The scaled Drude weight ZD for spinless fermions with 1 /Z-terms, obtained by evaluating 
the r.h.s. of (20) in d = 3 (i.e. Z = 6) with the actual three-dimensional DOS N (s) and 
N c fi(s), is shown in Figs. 1,2 as a function of interaction U and doping S, respectively. Also 
included in Fig. 1 is the total weight Z J dua(co) at 6 = which is essentially equal to the 
kinetic energy Z dua(uj) = —2ire 2 {H kin ) (/-sum rule) [33,34]. The difference between 
the two curves for 6 = represents the weight outside the <5-peak, i.e. duja^^uj). 
Both ZD and the total weight decrease quadratically with U at small U (Fig. 1). A large 
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U, however, ZD oc (hit/) 1 which implies a very slow decrease. The latter behavior is a 
consequence of the suppression of long-range order: in the symmetry-unbroken phase J] Fock 
adds to the hopping such that extended states are favored even at large U. 
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FIG. 1. Scaled Drude weight ZD vs. interaction U for spinless fermions with 1/Z-terms at 
T = for different fillings (D in units of e 2 t* j (fia d ~ 2 ) , where a is the lattice constant; U in units 
of t*). Solid curve: total weight § dwo{oj) (same units as D) for 6 = 0. 

While the decrease of the Drude weight with increasing repulsion is to be expected, its 
^-dependence - namely the global suppression of ZD upon doping - needs explanation. The 
latter effect is seen most clearly in Fig. 2 where ZD is found to decrease monotonically 
as a function of 6. This is surprising since close to half filling (small 6) one would expect 
doping to improve the mobility of the particles, as it is observed in the case of the Hubbard 
model [2,4,6]. Here the essential difference between the interaction in the Hubbard model 
and in the spinless fermion model comes into play. Let us consider a half-filled band at 
not too small U. In the Hubbard model the lattice sites are then singly occupied. The 
particles are thereby essentially localized, i.e. the Drude weight D vanishes. This effect is 
independent of whether the spins are antiferromagnetically ordered or not. Spin ordering, 
through which the energy is lowered by a (small) amount ~ t 2 /U, is here an additional 
effect with rather little influence on D, since D is already essentially zero. By contrast, for 
spinless fermions such a localization effect, and hence ZD ~ at 6 = 0, can only occur 
in the presence of long-range charge-order (checker board structure on the two sublattices) 
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since then the nearest-neighbor interaction is effectively avoided. In this case doping would 
indeed improve the mobility and thus increase ZD. In the homogeneous phase (which we 
consider here precisely because we wish to ignore the obvious effects of long-range order; 
see the discussion at the end of Section I) this localization effect is absent and hence ZD 
always decreases as a function of 6. At 6 = 0.29 there is a weak cusp due to the van Hove 
singularity of the DOS in d = 3 at one third of the band width. Note that in the empty band 
limit (6 — ► 0.5) ZD vanishes with the same slope for all U since the effect of the interaction 
decreases with decreasing particle density. 
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FIG. 2. Scale Drude weight ZD vs. doping 6 for spinless fermions with 1/Z-terms at T = for 
different U values; same units as in Fig. 1. 

B. T > 0: study of the dc-conductivity 

At low temperatures the dc-conductivity for the Hubbard model is found from (7), (13), 
(14) as 

z „ } --lfy^JMy (2I) 

where N^(uj) = — 7r~ 1 ImS JJ (a;). 

In analogy, for spinless fermions with 1/Z-terms one obtains 
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2 2 

ry spinless /rp\ ^ T 

/j<J dc K 1 ) — 

7T 



N c (co) , df(u 



(22) 



Here we encounter a peculiarity due to the fact that in the case of spinless fermions the 
limits Z — » oo and to — » do not commute: although in the limit Z — > oo the dynamic 
conductivity Z<j^ n,e "(w,T) - 0(1) for all to ^ 0, and also J? dua spinless (u, T) ~ 0(1), 
t/ie dc- conductivity behaves differently since N^(co) ~ 0(1/ Z). Equ. (22) then implies that 
^^"''"(T) - 0(Z) for Z ^ oo, i.e. the resistivity obeys p »^»»(T)/Z ~ (1/Z), while 



Hubbard 



(T)/Z ~ 0(1). 



The results for p(T)/Z in the case of spinless fermions with 1/Z-terms at half filling, 
obtained by evaluating the inverse of the r.h.s. of (22) in d = 3 (Z = 6), are shown in 
Figs. 3 and 4. They are compared with the results of Pruschke, Cox and Jarrell [10] for 
the half-filled Hubbard model in d = oo, obtained in the non-crossing approximation. In 
both cases U = 4.243. (We note that for the Hubbard model this value of U is just below 
the value where the Mott-Hubbard transition occurs [25,26,10]; the system is therefore quite 
close to the transition). 
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FIG. 3. Scaled resistivity pjZ vs. temperature T at U = 4.243t* and 6 = (p in units of 
ha d ~ 2 /e 2 ;T in units of t*). Main figure: spinless fermions with 1/Z-terms in d = 3. Inset: 
Hubbard model in d = oo, evaluated in the non-crossing approximation [8]; see text. 
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FIG. 4. Same results as in Fig. 3 on double logarithmic scale; same units as in Fig. 3. The 
T 2 -behavior at low T is clearly seen. 

At low temperatures, T ~ 0.1, both models show Fermi liquid behavior, i.e. p(T)/Z = 
AT 2 , as is clearly seen on a double-logarithmic plot (Fig. 4) [35]. The typical T 2 -dependence 
originates in the low-frequency behavior of N^(ui) = —ir~ 1 Im'E R (uj) =: 7V S (0)(T 2 + bco 2 ) in 
(21), (22) where b is some constant. Since (—df/du) is a strongly peaked function of ui at 
the Fermi level, with width ~ T, one obtains p(T)/Z = (Za^y 1 oc iV s (0)T 2 i.e. A oc %(0). 
For spinless fermions with 1/Z-terms one has iV s (0) oc U 2 /Z for all U. For the Hubbard 
model iV s (0) can become much larger (this is so even at small U, where iV s (0) oc U 2 in this 
case), in particular close to the metal-insulator transition. This explains why A is much 
larger for the Hubbard model (see Figs. 3,4). 

The existence of a maximum in the resistivity for the Hubbard model is in striking 
contrast to the monotonic increase of p(T)/Z in the case of spinless fermions with 1/Z- 
terms. This maximum is located almost exactly at the temperature where the system would 
undergo a phase transition to the antiferromagnetic phase [11] if that ordering were not 
suppressed. Hence the increase in the resistivity can be attributed to the enhancement in 
the scattering of electrons by local spin fluctuations. Such an effect does not occur in the 
case of spinless fermions with 1/Z-terms: the interaction U jZ is by a factor \ jZ weaker and 
the charge fluctuations are averaged over the very large number of nearest neighbors, i.e. do 
not act coherently. 
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At high temperatures the resistivity increases linearly with temperature in both models 
(for the Hubbard model that part is not shown in Figs. 3,4). This is due to the fact that 
for T — » oo the functions N(u), N c (u>) and N^(u) become temperature independent, while 
(-df/du) -> 1/(4T), thus leading to Za dc oc 1/T, i.e. p(T)/Z oc T. 

IV. DISCUSSION 

We will now discuss the conditions under which the dc-conductivity can be finite at 
T > 0, and to what extent this is a generic property of correlated electron systems in 
dimensions d = 1, 2, 3. 

A. Umklapp scattering 

In the absence of phonons and impurities only Umklapp-scattering processes are able 
to degrade the total crystal momentum. Under such conditions a finite dc-conductivity at 
T > can therefore only occur in a lattice system. In the d = oo dynamic mean-field theory 
for the Hubbard model and spinless fermions with 1/Z-corrections, the correlation problem 
reduces to a self-consistent single-site problem [22-24], where the importance of the lattice 
and the notion of reciprocal lattice vectors are no longer directly visible. Nevertheless it 
must be borne in mind that for quantum mechanical particles the d — ► oo limit can only be 
formulated on a lattice (the lattice constant provides a natural short wavelength cut-off [15]). 
Hence the results obtained for correlated elecrons in d = oo, in particular those obtained 
here for the Drude weight D and a^c, require the existence of a lattice. 

In fact, in d = oo all scattering processes involve Umklapp processes. This is a direct 
consequence of the topology of the Fermi surface of a lattice system which in the limit d — >• oo 
is very close to a sphere from which major parts are chopped off at the Brillouin-zone (BZ) 
boundaries [36]. In fact, the Fermi body covers the fraction n/g of the (d — l)-dimensional 
boundary surface of the hypercubic BZ, where n is the particle density and g denotes the 
degeneracy of the particle state, i.e. g = 1 and 2 for spinless fermions and the Hubbard 
model, respectively. Since in d = oo the Fermi body intersects with the BZ boundary for 
all fillings n it is possible to have Umklapp processes with arbitrarily small q '-vectors. There 
is then no minimum energy A required for Umklapp scattering and hence no exponential 
suppression exp(— A/T) of such processes at low temperature. Thus their contribution is 
always important. 

For each dimension d of a hypercubic lattice one may define a critical density n c ^ such 
that for n > n[ d ^ the Fermi body intesects with the BZ boundaries. While n Cj00 = (note, 
however, that the limits n — » and d — > oo do not commute), one has < n c ^ < 1 for 
any finite dimension d with 3 < d < oo, i.e. in general the critical density increases with 
decreasing d. In particular, in d = 2 one has n Cj2 = 1, in which case the Fermi body touches 



12 



the BZ at discrete points, while in d = 1 an intersection does not occur for any filling. 
Hence Umklapp processes become more effective the higher the dimension of the lattice is. 
Consequently the d = oo mean-field theory overestimates the effect of Umklapp processes. 
A finite dc-conductivity is nevertheless a natural feature of lattice systems in sufficiently 
high dimensions d. 

B. Vertex corrections 

In the limit d = oo the dynamic conductivity a(u, T) is strictly given by the dressed 
bubble-diagram since vertex corrections are suppressed [29]. The dc-conductivity <r dc (T) is 
thus determined by the single- particle life time. In general the life-time (or rather its inverse, 
the scattering rate) receives contributions from scattering processes of quasi-particles involv- 
ing Umklapp scattering as well as momentum-conserving "normal" scattering. In d = oo 
essentially all scattering events involve Umklapp processes. Genuine two-particle interference 
effects are described by vertex corrections. They enter as soon as 1/Z-corrections are taken 
into account and lead, for example, to the renormalization factor 7 = 1 + O(^) > 1, (18), 
in the case of spinless fermions. By including 1 jZ corrections, i.e. vertex corrections, one 
departs from the d = 00 limit. Thereby the importance of the current-degrading Umklapp 
processes is reduced. Hence we expect 1 /Z-corrections to increase the conductivity. 

It is still an unresolved question whether or not vertex corrections are able to restore 
momentum conservation to such an extent that cr^ c (T) can be infinite at T > 0, i. e. whether 
<j(u,T) has a contribution D(T)6(u) even for T > 0. If at all, this would be a property 
of systems where Umklapp processes are relatively unimportant, such as Zoty-dimensional 
system away from half- filling [9]. Recently Castella et al. [13] conjectured that a divergent 
dc-conductivity at T > is a generic property of integrable system such as the Luttinger 
model, interacting bosons and the U = 00 Hubbard model in d = 1. In this case our finding 
of a finite dc-conductivity at T > would be generic in any dimension since most interaction 
models appear to be non-integrable. 

It will be very interesting to investigate further whether, and if so under what conditions, 
the dc-conductivity of correlated lattice electrons may diverge at T > even in dimensions 
d = 2 or 3. 

We are grateful to H. Fukuyama for generating our interest in the problem, to G. Czy- 
choll, W. Metzner, E. Miiller-Hartmann for valuable discussions, and to W. Hanke, 
P. Horsch, S. Maekawa and Th. Pruschke for useful information. This work was supported 
in part by the SFB 341 of the Deutsche Forschungsgemeinschaft. 
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